Rabi Oscillations and Few-Level Approximations 
in Time-Dependent Density Functional Theory 
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The resonant interaction of laser light with atoms is analyzed from the time-dependent density 
functional theory perspective using a model Helium atom which can be solved exactly. It is found 
that in exact-exchange approximation the time-dependent dipole shows Rabi-type oscillations of its 
amplitude. However, the time-dependent density itself is not well described. These seemingly con- 
tradictory findings are analyzed. The Rabi-type oscillations are found to be essentially of classical 
origin. The incompatibility of time-dependent density functional theory with few-level approxima- 
tions for the description of resonant dynamics is discussed. 

PACS numbers: 31.15.ee, 31.70.Hq 



Linear response time-dependent density functional the- 
ory is widely and successfully applied to calculate absorp- 
tion spectra of atoms, molecules, and solids [![. The pre- 
dicted transitions between the groundstate and singly ex- 
cited states are often remarkably accurate Q, even with 
simple approximations to the exchange-correlation po- 
tential. In experiments, a resonant interaction between 
a laser field and atoms or molecules is routinely used 
to, e.g., prepare the system in an excited state. The 
time it takes to transfer the groundstate to an excited 
state is half a Rabi period (7r-pulsc) and can (in sim- 
ple cases) be calculated using a two-level approximation 
(TLA, see any text book on quantum optics, e.g., 0]). 
Time-dependent density functional theory (TDDFT) be- 
yond linear response is, in principle, capable of describing 
the entire dynamics of the electron density n(r, t) exactly 
if the exact exchange-correlation potential were known 
With the increasing interest in real-time quantum 
dynamics of matter exposed to laser light TDDFT be- 
yond linear response attracts more and more attention 
(see [l| and references therein). 

In this Letter we analyze the resonant interaction 
of laser light with atoms from the TDDFT perspec- 
tive. To that end we employ a numerically soluble one- 
dimensional helium atom as a benchmark model for the 
corresponding time-dependent Kohn-Sham (TDKS) cal- 
culations. We show that the TDKS dipole indeed dis- 
plays Rabi-type oscillations which, however, are of clas- 
sical origin, and that the density itself is not properly de- 
scribed in the exact exchange-only approximation. The 
incompatibility of a TLA with the TDKS equation in the 
case of resonant interaction is also discussed. 

Let us first introduce the widely used one-dimensional 
model helium [|| in which both electrons move along the 
laser polarization direction only. Softening the Coulomb 
interaction, l/|r| — > l/\/l + x 2 [atomic units (a.u.) are 
used], we obtain for the Hamiltonian H(t) = T+V C xt{t) + 
V ee with V ext (t) = -2/Vl + x 2 -2/\/l + x' i +E{t)[x + 
x'], T = -(l/2)[fl£ + dl,\ and V ee = 1/y/l + {x - x'f. 
Here, the coupling to the laser field is described by 
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FIG. 1: (Color online) Linear response of the He model atom 
as obtained from the full TDSE (black, labelled 'TDSE'), from 
the x-only TDKS (red, labelled 'TDKS'), and the frozen x- 
only KS potential (green, labelled 'frozen'). The vertical ar- 
row indicates the first ionization threshold. The insert shows 
a close-up of the transitions to singly excited states. 



the electric field in dipole approximation E(t), T is 
the kinetic energy, and V ee the electron-electron inter- 
action. We start from the spin-singlet groundstate. 
Since there is no spin-dependent external potential the 
state will remain a spin-singlet state, i.e., (xx'\^(t)) = 
(1/V2)ip(xx't)[\ U) - I IT)] with ip(xx't) being symmet- 
ric under the exchange of x and x' . The time-evolution 
is governed by the time-dependent Schrodingcr equation 
(TDSE) id t if(xx't) = H{t)ip(xx't). The groundstate en- 
ergy in this system is £q = —2.238, the first excited spin- 
singlet state is at £± = —1.705. 

The corresponding TDKS equation reads idtf(xt) = 
[-(1/2)^2 + v KS (xt)](p{xt) with 



E(t)x 



\<p{x't)\ 2 dx' 
^l + (x~x') 2 



(1) 



in exact exchange-only approximation, the last term be- 
ing the Hartree-exchange potential vu x . In this ap- 
proximation correlation effects are neglected, and the 
TDKS equation equals the time-dependent Hartree-Fock 
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(TDHF) equation. The TDHF wave function is a Slater- 
determinant, {xx'\<f>(t)) = ip(x't)tp(xt)[\ U) - | 4T)]/\/2, 
and the TDKS density is simply given by n(xt) = 
2\cp(xt)\ 2 . 

The linear response spectra are calculated according to 
Ref. @. The result is shown in Fig. Q] The strongest peak 
is associated with the transition between the groundstate 
and the first excited singlet-state at co = E^—Eq = 0.533, 
followed by transitions to higher excited states and the 
first continuum (indicated by the vertical arrow). The 
TDSE-spectrum also shows transitions to doubly excited 
states and the corresponding continua. Such transitions 
are absent in linear response TDDFT employing simple, 
adiabatic exchange-correlation potentials 0] ■ 

The linear response spectrum for the "frozen" 
Kohn-Sham (KS) groundstate potential Vks^) = 
-2/Vl + x 2 + J in {x') dz'VCV 1 + {x~ x 1 ) 2 ) (com- 
monly called "bare" KS response) is also included in 
Fig. [TJ showing the transitions to excited states in v^ s (x) 
and illustrating that VKs{xt) — v^(x) shifts the peaks 
closer to the correct positions. 

Let us now consider a monochromatic laser beam of 

T I 

resonant frequency to = £{ — £q and electric field ampli- 
tude E = io A. Assuming that the TLA and rotating wave 
approximation are valid, we expect Rabi oscillations of 
frequency f2 = -E/iio to occur, where 

Mio = (^oK 2 -' + ^Ol^i) = 2 J J Ax dxipQ(xx')xipi{xx) 

is the transition dipole matrix element. Its numerical 
value is 1.1. The dipole then evolves according 

d(t) = fiio sin cot shiQt. (2) 

Figure [2] shows (x)(t) = d(t)/2 as it results from the 
TDSE, TDKS, and frozen KS calculations. A laser field 
of vector potential amplitude A = 0.0125 was ramped 
up over two laser cycles and then held constant. The 
laser frequency was tuned to the resonance to = £{ — £q, 
i.e., lo = 0.533 for the TDSE, lo = 0.549 for the TDKS, 
and lo = 0.492 for the frozen KS calculation (all in- 
ferred from Fig. [T|). The TDSE result shown in panel 
(a) displays Rabi oscillations of the envelope of frequency 
tt = lw/iio = 0.0075, as expected. At t = ir/Q, ~ 420 
the excited state is maximally populated and the enve- 
lope of the excursion is close to zero. At t = 2tt/H ~ 840 
the system is mostly in the groundstate again. A closer 
inspection of the TDSE result shows that because of ion- 
ization and transitions to other states the population of 
the first excited state after half the Rabi period is only 
0.975 instead of unity. The population of the ground- 
state after a full Rabi cycle is 0.96. In the frozen KS 
calculation [panel (c)] ionization and the population of 
other excited states are more pronounced. As a conse- 
quence, the excursion envelope does not go to zero at 
t = 7r/Sl and the excursion amplitude is overestimated. 
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FIG. 2: (Color online) Expectation values (x) vs time for a 
resonant excitation with A — 0.0125 and (a) lo — 0.533 for 
the TDSE, (b) lo = 0.549 for the TDKS, and (c) lo = 0.492 
for the frozen KS calculation. 
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FIG. 3: (Color online) Exact groundstate density no(x)/2 = 
J dx' ^(xx')! 2 (black, labelled 'ground') and exact excited 
state density m(x)/2 = J dx' \ipi(xx')\ 2 (red, labelled 'ex- 
cited'). The two inserts show contour plots of |?/>o(xx')| 2 (left) 
and IV'i(xx')! 2 (right). The signs of ipi(xx') are indicated in 
the right insert. 

The TDKS calculation in panel (b) shows oscillations of 
the right amplitude. The Rabi period one infers from 
the envelope-oscillations is remarkably close to the ex- 
act result in (a). For non-resonant driving the amplitude 
oscillations are absent, as they should. 

Figure ^> suggests that Rabi oscillations are well de- 
scribed within the TDKS system. If this were true the 
TDKS density should oscillate between the groundstate 
density and a density similar to the exact first excited 
state density shown in Fig. [3] Unfortunately, this is not 
the case. Examining the TDKS density at time t ~ 350 
reveals that it does not assume the shape of the exact 
excited state density of Fig. [3] but rather resembles the 
initial density again! Hence, despite an erroneous TDKS 
density we observe the Rabi-like oscillations of Fig. [2Jd 
in its first moment, i.e., in the TDKS dipole. As for our 
two-electron system the exact KS orbital corresponding 
to the excited state density is given by ipi (x) = y/ni(x)/2 
and thus, according to Fig. [31 has no nodes, the exact 
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KS orbital representing this excited state density must 
be the groundstate of a KS potential wks["i]- Hence, 
even the exact KS potential will not lead to a population 
transfer to an excited KS state but will "guide" the den- 
sity towards the stationary, "new" groundstate density ri\ 
during a 7r-pulse. The exact exchange-only approxima- 
tion used in ^ above does not do this. Hence, correla- 
tion is needed to describe the density dynamics properly. 
Whether memory effects Q are important in this con- 
text will be investigated in a forthcoming paper. Note 
that the exact KS potential may be simply constructed 
by inversion for the case of the two-electron spin-singlet 
system studied here Q. However, in order to actually 
identify memory effects more involved methods are re- 
quired, e.g., the approach proposed in Ref. [loj . 

Let us now investigate the origin of the Rabi-like oscil- 
lations in the TDKS dipole of Fig. [2] Since it is not due to 
the density-dynamics corresponding to population trans- 
fer there must be another explanation. In order to show 
that the oscillations are classical in origin let the motion 
of the center of mass of the density being approximated 
by the dynamics of a point particle in an effective an- 
harmonic potential v(x) = (ui 2 /2)x 2 + (a/3)x 3 + (/3/4)x 4 
with a and f3 some constants. The external driver is 
of the form E(t) = Ecos[(ui + e)t], i.e., e is the detun- 
ing with respect to u>, i.e., the frequency characterizing 
the harmonic region of the potential around the origin. 
The squared excursion amplitude B 2 (e) = [maxi(t)] 2 



fulfills a cubic equation [111 ]. As soon as a critical driver 



strength is reached more than one real solution for B(e) 
exists, and a discontinuity develops. Figure |4] shows that 
the amplitude of the TDKS dipole as a function of the 
laser frequency displays exactly this feature. The clas- 
sical Rabi-likc oscillations can intuitively be understood 
as follows: while being in the harmonic region of the 
potential the particle is resonantly driven and thus the 
excursion amplitude increases. However, with increasing 
excursion amplitude the particle inevitably senses the an- 
harmonicity of the potential. As a consequence the driver 
is not resonant anymore and the excursion amplitude de- 
creases. Hence the Rabi-like oscillations of the excur- 
sion amplitude seen in our TDKS results arc essentially 
of classical origin. Similar oscillations were observed in 
Josephson junctions 

Going back to our model helium where 
the TLA applied to the interacting system, 
|tt(t)} ~ a(i)exp(-i£ i)|*o) + &(*) exp(-i£i*)|*i>, 
accurately captures the resonant population trans- 
excited state, the density evolves in time 



fer to an 
according 



n(xt) 



\a(t)\ 2 n (x) + \b{t)\ 2 ni {x) 



2$.{a*(t)b(t)exp(-iojt)An(x)}. Here An(x) = 

(\l>o|7i(x)|\f r i) is real with h(x) = J^a H{ x )^<j( x ) 
the density operator in second quantization with spin 
degrees of freedom a. Can we use this accurate, 
interacting two-level density to construct a i>Hxc for 
a corresponding two-level KS scheme? In order to 




0.50 0.55 0.60 

Frequency (a.u.) 

FIG. 4: TDKS excursion amplitude B — maxi vs laser fre- 
quency for a driver with A = 0.0125. Discontinuity and 
asymmetric peak structure are characteristic of classical an- 
harmonic oscillations lllll. 



show that such an approach will fail we start from the 
fundamental equation for the derivation of the extended 
Runge-Gross proof (see [IH or Ch. 2 in in one 
dimension, 

d 2 n(xt) = d x [n(xt)d x v cxt (xt)] + q([n];xt), (3) 

with q([n];xt) = (*(t)\d 2 x f(x) + 0x%, (*)!*(*)>, 
the momentum-stress tensor T(x) = 
E CT {[ 9 x^( x )]<9xVv(a;)- jd 2 [4>t{x)ipa(x)]} and the inter- 
action term V^ c (x) = Y.a,a> I dx ' i>t( x )'$t> {x')[d x v cc {\x- 
x'\))ifja<(x')ija{x). With the analogue of ([3]) for 
the noninteracting system and its state |3>(i)) 
[leading to the same density n(xt)], one finds with 
wks(M; srf) = v C xt(xt) + whxc ( [n] ; xi) for the Hartree- 
cxchangc-correlation potential 

d x [n{xt)d x VK XC ([n]-,xt)} = q([n];xt) (4) 

-m)\d 2 x T(x)mt)). 

One could now use the interacting two-level density in 
(|3|) to construct ^xc([n];xt). However, even if one 
were able to determine a{t) and b(t) from the nonin- 
teracting system alone (so that one has not to solve 
the full interacting problem in the first place) such a 
TDKS treatment will lead to wrong predictions. Note 
that with ([3]) and v ex t( x t) — vq(x) + xE{t) the dipole 
acceleration of an N particle system reads d(t) = 
—NE(t) — J dx n(xt)d x vo(x). Here the term depending 
on d x vn xc ([n\; xt) vanishes in accordance with the zero- 
force theorem [l], Q|. With n Q (x), ni(x) symmetric and 
Vq(x) even, J dx n a (x)d x Vo(x) = j dx ni(x)d x v Q (x) = 
results, and thus 



d(t) + cd(t) = —NE(t), 



(5) 



where c = J dx (d x vo(x))An(x)/ J dxxAn(x). Equation 
([5]) describes a driven harmonic oscillator which [for ini- 
tial conditions d(0) = d(0) = 0] does not exhibit oscil- 
lations of the excursion amplitude as a function of the 
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driver amplitude, i.e., no Rabi-like oscillations. Hence, 
despite using an accurate density as input not even the 
dipole is reproduced within such a two- level TDDFT. 
What went wrong? It turns out that the introduction of 
a TLA into leads to inconsistencies between Hilbert 
spaces and Hamiltonians. In order to illustrate this fact 
it is sufficient to consider the Heisenberg equation of mo- 
tion for some general, time-independent operator O, 



d t On(t) = i[Hn(t),6 H (t)] 

= i(u- l {t)[H{t),6]U(t) 



(G) 



with 6 H (» = U^itfOUit) and U(t) 
T(exp(-i/ *di'F(t'))) the time evolution opera- 
tor. Making use of ([3]) and introducing a TLA amounts 
to calculating the commutators in expressions like ([6|) 
in the full Hilbert space while the time evolution is 
performed in a reduced Hilbert space, i.e., in our case in 
a two-level subspace, 



i (U^(t)[H(t),0]U 2 (t) 



(7) 



i{Uz\t)i 2 (H{t)i6-diH{tj) i 2 u 2 (t)} 



with i = E^ol'Wfcl, la = £Ul*kX**l and 
U 2 (t) = T (cxp(-i J * dt' i 2 H(t')i 2 )y Instead, for a 
consistent TLA 

d t (u^(t)OU 2 (t)) (8) 

= i {u 2 Ht)U [H{t)i 2 d - 6i 2 H{t)) i 2 i> 2 (<)} 

holds, with a different commutator leading to a differ- 
ent equation of motion. In particular, the equation for 
dfn(xt) is different from ([3]) if the Hamiltonian is re- 
stricted to a two- level subspace [15| . Obviously, our anal- 
ysis not only applies to a TLA but to any finite-level 
approximation. 

So far we tried to construct uh xc by applying a TLA 
to Eq. On the other hand, by inversion [§] we 

can determine the KS potential «ks generating exactly a 
given two-level density n(xt). The associated dipole ac- 
celeration then reads d(t) = — J dxn(xt)d x VKs([n];xt). 
Subtracting from the KS potential the physical ex- 
ternal potential defines «hxc, and we obtain d(t) + 
J dxn(xt)d x vo(x) + f dx n(xt)d x vn xc ([n}; xt) = d(t) + 
cd(t) + J dx n(xt)d x V}i xc ([n]; xt) = —NE(t). By con- 
struction, the two-level density leads to the correct 
two-level dipole acceleration. This can only be possi- 
ble if the term depending on vhxc contributes. Oth- 
erwise the same problem as with Eq. ([5]) discussed 
above arises. However, a nonvanishing contribution from 
/ dx n{xt)d x v^ c {[n]\ xt) is only possible if i?hxc does not 
describe internal forces only and thus violates the zero- 
force theorem [l|, [lj. As a consequence, the external po- 
tentials of the interacting and the noninteracting system 



cannot be kept equal. In fact, there is no local external 
potential v ex t(x) that supports just two levels. Neverthe- 
less, a uhxc, which violates the zero-force theorem, may 
be acceptable as an approximation. 

In conclusion, we investigated resonant dynamics in 
the exact- exchange approximation for the simple but 
numerically exactly solvable case of a one-dimensional 
model helium atom. Although the dipole shows Rabi- 
type oscillations the density-dynamics of the popula- 
tion transfer process is not properly described. As 
a consequence, the dipole spectra calculated using ex- 
act exchange-only TDDFT may be sufficiently accurate 
while the real-time dynamics of the density is erroneous. 
The incompatibility of few-level approximations with 
TDDFT to describe resonant density-dynamics was an- 
alyzed. Since there is an increasing interest in strongly 
driven real-time quantum dynamics of matter the devel- 
opment of exchange-correlation potentials capable of de- 
scribing resonant charge transfer is particularly impor- 
tant and will be the subject of future work. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft. 
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